G^ 



o 



0^ 



gr-qc/9712043 

DAMTP Relativity Group Preprint: DAMTP-R97/56 

Fermilab Preprint: Pub-97/364-A 



Persistence Amplitudes from Numerical Quantum 
g;": Gravity 



p. D'Eath, A. Sornborger, 



^ ; DAMTP, NASA/Fermilab Astrophysics Group, 

Cambridge University, Fermi National Accelerator Laboratory, 



3 Silver Street, Box 500, 



> ■ Cambridge, CB3 9EW, Batavia, IL 60510-0500 

^ ■ UK USA 

o 

<N ■ October 27, 1997 

O . Abstract 

^ [ The Euclidean quantum amplitude to go between data specified on an initial 

p|J' and a final hypersurface may be approximated by the tree amplitude 

k> ; exp{-Idassical/fi), (1) 

C^ '_ where I classical is the Euclidean action of the classical solution joining the initial 

and final data. In certain cases the tree amplitude is exact. We study Idassicah s-nd 
hence the quantum amplitude, in the case of a spherically symmetric Riemannian 
gravitational field coupled to a spherically symmetric scalar field. The classical 
scalar field obeys an elliptic equation, which we solve using relaxation techniques, 
in conjunction with the field equations giving the gravitational field. An example 
of the transition from linearity to non-linearity is presented and power law behavior 
of the action is demonstrated. 

PACS numbers 04.20.-q, 04.20.Jb, 04.60. +n 

1 Introduction 

Much can be learnt in quantum gravity by considering the amplitude to go from an initial 
to a final configuration of a gravitational and a massless scalar field 0. To complete 





the specification of tlie amplitude, one needs to give tlie time separation between tlie 
two spacelike surfaces depicted, together with other data at spatial infinity [1]. If, for 
example, the time separation is taken to be large, and the bounding data are taken to 
be weak, then one has a scattering configuration. On the the other hand, the formalism 
is equally well adapted to studying strong-field amplitudes, as we shall see in this paper, 
with the help of numerical analysis. One such example concerns the late-time evolution 
of a radiating black hole. Further examples are provided in quantum cosmology [2-4], 
where one again needs the strong-field description. 

Formally (if the time separation is taken to be Lorentzian) the amplitude is given by 
a Lorentzian path integral 

amplitude = f e^'^^^\ (2) 

Here the fields g^u and (p in the integration must agree with the boundary data, the 
intrinsic three-dimensional metric hij and the scalar field 0. Quite apart from convergence 
questions, there is a more elementary difficulty with the Lorentzian path integral. One 
might expect the path integral (0) to admit a semi-classical expansion 

amplitude ~ (A + hAi + ■ ■ ■)exp{iSc/h). (3) 

Here Sc denotes the classical action (if such exists) of a classical solution joining the 
initial to the final data. It is a functional of the boundary data and data at infinity. 
Similarly the one-loop term A, two-loop term Ai and so forth depend on the boundary 
data. Now consider the special case in which the classical geometry is (nearly) flat and 
one has a massless scalar field between an initial surface ti = where = and a final 
surface at time ^2 where is non-zero [4]. The classical solution is 

0(t,x) = const, f d^k^^^i^e'''-''stn(\k\t). (4) 

For general final data, the 'solution' is singular because of the poles in the denomina- 
tor. This is an example of a well-known phenomenon, that the boundary- value problem 
for a hyperbolic system of equations is not well posed [5]. The expression (|) does, how- 
ever, make it clear that the boundary- value problem is well-posed if the time-separation 
is rotated by a small amount ie into the complex plane[6]. 

The most natural arena for studying boundary- value problems is within Riemannian 
geometry, where the Lorentzian time coordinate t is replaced by t = —ir, with a positive- 
definite four-metric g^u and massless scalar field [7]. One might conjecture that the 
classical boundary- value problem is well-posed for moderate-sized boundary data; indeed, 
the numerical methods to be described in Sec. 2 will test this conjecture. The conjecture 
has been shown analytically to hold for weak (but non-linear) gravitational fields in the 
absence of a scalar field [8]. The full question of the existence of the coupled nonlinear 
elliptic Einstein-scalar field equations subject to boundary data is a major unresolved 
problem in general relativity and partial differential equations. As suggested in the 
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previous paragraph, one would expect to obtain a Lorentzian amplitude by starting 
with a Euclidean time interval at infinity, and then rotating the Euclidean time interval 
towards a Lorentzian time interval. The first treatment of the black-hole evaporation 
problem by Euclidean methods was by Hartle and Hawking [9]. The original Euclidean 
approach to quantum cosmology was also due to Hartle and Hawking [3]. 

For a Riemannian classical geometry, the quantum amplitude will admit the semi- 
classical expansion 

amplitude ~ (A + hAi + ■ ■ ■)exp{—Ic/h). (5) 

Here Ic is the classical action for gravity coupled to a massless scalar, and A,Ai,--- 
are one-, two-, and higher-loop terms, dependent on the boundary data. Even at the 
one-loop level, for pure gravity, the Riemannian amplitude is known to be divergent, 
with divergence linear in the surface invariants [4,10] 

/2 = jd^xh^l'^itrKf, 

Jg = jd^xh^/'^Ki^K'^itrK), ^^^ 

h = Jd^xh^/^K'jKiK^. 

Here Rij is the Ricci tensor formed from the three-dimensional metric hij, and Kij is 
the second fundamental form [11]. (Without boundaries, pure gravity is divergent at two 
loops and beyond [12].) When a massless scalar field is added, there are again one-loop 
surface divergences. 

All these divergences may be cured by including supersymmetry, providing that one 
continues to work with surface boundary data given by hij in the purely gravitational 
case, and with hij,(f) in the gravity- scalar case. From this point of view quantum cor- 
rections are 'soft' in supersymmetric theories. In fact, provided one sets the appropriate 
supersymmetric partners to zero on the boundaries, it turns out that the amplitude just 
has the tree form exp{—Ic/h), where Ic is the classical action of the original bosonic 
fields. One may use directly the supergravity calculations of [4]. Note that there was a 
mistake in the paper [13] which included an earlier discussion of the ideas of [4]. This has 
been corrected in [4] by the inclusion of the auxiliary fields of supergravity [14] , which are 
of course needed for a complete treatment of supergravity. Another way of motivating 
this calculation is to understand the action as being derived as an effective action from 
low-energy superstring theory, in which case we are calculating tree-level amplitudes. 

The local boundary conditions used here have been stressed in part because they 
appear naturally in specifying the path integral. However it is also important to under- 
stand that scattering boundary conditions, by contrast, may have very different proper- 
ties. One obtains scattering boundary conditions by pushing the two spacelike surfaces 



off to infinity. Then one takes an infinite product of zero- or one- occupancy harmonic 
oscillator states. This implies that there is a non-trivial path integral involved in mov- 
ing between local and scattering boundary conditions, so that it is not surprising that 
quantum amplitudes with the two types of boundary conditions should be very differ- 
ent. Indeed [15], with scattering boundary conditions the supersymmetric gravity-scalar 
model is divergent at one loop. 

Since the quantum amplitude for the (supersymmetrized) gravity- scalar model is ex- 
actly semi-classical, one strategy for computing such amplitudes is clear: set up boundary 
data and solve the coupled field equations numerically. Include the classical Euclidean 
action I^ in the numerical calculation, and finally compute the amplitude exp{—Ic/h). 
The equations for the fields and an expression for the action are presented in Sec. 2. In 
Sec. 3, we discuss the numerical methods used to find solutions. In Sec. 4, we present 
the results obtained computationally. And, in Sec. 5, we discuss the results and present 
our conclusions. 

2 Spherical Riemannian Scalar- Gravity Solutions 

The Riemannian metric is taken for simplicity in the spherically symmetric form 

Qf,^ = diag{e'', e"-, r^, r^sin^9), (7) 

in coordinates {t,r,6,(j)), where t is a Riemannian time coordinate and 

a = a{t,r), b = b{t,r). (8) 

Notice that for real values of a and b in the metric, the metric is positive definite. 
The scalar field is also taken to be spherically symmetric: 

= 0(t,r). (9) 

And the Euclidean action has the form 

/ = -j^Jd^xg^R 

+ \^d'xg'^{V(^f 

(10) 

+ boundary contributions. 



Here g = det{g^y) and R is the Ricci scalar of g^j^y. The boundary contributions 
will later be important and are discussed in Eqs.(2.19-22). The field equations may be 
written in the form 



d,{gh^''<P,.) = ^^^^ 

Here, we have set G = ^ = c = 1. 

The Lorentzian version of this problem, with a spherically symmetric scalar field 
in a Lorentzian spherically symmetric gravitational field, has been studied extensively, 
leading to many interesting results [16-22]. 

Explicitly, the Euclidean gravitational field equations are 

Rtt = ie-"+V6' - |a - id^ 
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Rrr = ia'b' + ^ - ^e-'+^'d - {e-'+'^a 
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The field equations can be simplified, to give (starting with the scalar field equation): 

+ e'-y" + ^(d - 6)0 + ^(1 + e'^)0' = 0, (13) 
2 r 

a' = -47rr(e«-V' - 0'") + il^l^ (14) 

b' = -47rr(e'^-V' - 0'') - ii^l^ (15) 

r 

d = 87rr00' (16) 

a + e'-"6" + ^(d - b)d - e^-'^ (^~^") (5^ + ^^ = Sn{^^ + e'-"0''). (17) 



One can see the general structure of the coupled Riemannian Einstein-scalar equations 
from this system. The scalar field obeys an elliptic equation (|13]) which is determined 
if the gravitational background is known. Conversely, Eqs. (|14| - |1^) determine a and 
6 if is known. An effective nonlinearity is present in the elliptic Eq. (|13D due to the 
presence of a and h which involve 0. Notice that Eq. (|I^ is equivalent to Eq. (|T3|) given 
the constraints, i.e. they are not independent equations. 

As the simplest example of a boundary-value problem, one might wish to solve the 
field equations inside a rectangular boundary, where r runs from the axis r = to an 
outer boundary r = Vmax, and t runs from an initial value U to a final value tj. One 
normally expects that and the intrinsic 3-metric hij would be specified on the boundary 
[4], in which case one would have a = ai fixed on the initial surface, a^ on the final surface, 
and b = bf on the outer surface. Note also that, for regularity of the initial 3-metric on 
the initial surface, one should have a -^ as r — > on the initial surface. Now consider 
Eq. (|TB|). By regularity on the axis, 0' = there. Hence d = along the axis. Hence 

a(t,r = 0) = 0. (18) 

This provides another boundary condition. 

However, we are unable to solve this version of the boundary- value problem, as the 
data on the boundaries are overdetermined, we cannot simultaneously fix and a and 
b. The correct approach is to fix on the boundaries, as well as fixing b on the outer 
boundary (thus fixing the gauge freedom for b), then solve Eq. ([13|) by relaxation [23] 
and then to compute a and b by integrating the constraint equations (|1^, 0) iteratively, 
while checking that the full set (13 - 17) holds after the iteration. There is also one 
leftover gauge degree of freedom, the value of b, which we choose to fix at r = r^ax by 
fixing b at Vmax- 

The classical action for our system of equations resides on the boundary. It is 

Ic = / d^xhkrK + - f Sxh^cj)^. (19) 

Svr J 2 7 on 

The integral is taken over all bounding surfaces, with intrinsic metric hij and h = det{hij). 
Here Kij is the Euclidean second fundamental form [4], taken with respect to the outward 
normal, and trK = K'^Kij. For the upper surface t = tj (say), one finds 

If = --f drr^e^"-^)/^^ + 2tt f drr^e^^-^^/^cj)^, (20) 

with a corresponding expression for Jj in which the minus and plus signs are reversed. 
There is also a contribution from the outer boundary at r = rj, which is near r = oo in 
our case. It is [4] 

/oo = MT, (21) 

where M is the mass of the space-time and T is the proper Euclidean time-interval 
between the initial and final surfaces. Hence the total classical Euclidean action is 



total = -\ Udrr^e'---'^'^a + 271 j drr^^'^' 






3 Numerical Methods 

We use the standard artificial time method for relaxing the elhptic equation. This 
relaxation method uses the artifice of introducing a diffusive term to the elliptic equation. 
The diffusive term gradually relaxes to zero and one is left with a solution to the elliptic 
operator of interest. 

With this relaxation method, we solve the equation 

9^0 + L(0, a, 6, d, 6) = (23) 

where r is an artificial time coordinate, and L is the nonlinear elliptic operator for 
which we desire a solution. We solve the equation by discretising in the artificial time 
and also imaginary time and space, then integrating the equation in time. For a linear 
elliptic operator discretised in this manner, there is a constraint on the allowed size of 
the integration timestep. For instance, for solving the heat equation, the constraint is 
At < (Ax)^/2. For a nonlinear elliptic operator, the timestep must typically be kept 
substantially smaller, such that the integration is effectively linear. In our integrations, 
for a simulation volume of 0.8, we took At = 0.00005, typically. The small value of 
timestep needed in our case is a measure of the non-linearity of our problem. 

After one integration of Eq. ( PB| ) where L is given by (|r^), we use Eq. ( ]T^ ) to 
integrate a radially outwards from the axis where a = 0. For a radius of one gridpoint 
out from the axis we are forced to approximate Eq. (|14D in order to avoid problems from 
divergences due to small values of r. At large r, the solution for a will be dominated 
by the non-(/) part of Eq. (|14|). The solution to this will be a Schwarzschild solution, so 
that one will have a ~ 2M/r as r ~ cxo, where M is the mass of the space-time. One 
also integrates Eq. (|15D to derive values for b inwards from the outer boundary to the 
axis. Then one iterates this procedure until (f),a,b converge (if this occurs), checks that 
the remaining equations ( pJGJ p!7| ) hold, and evaluates the action. 

4 Results 

We have run a number of integrations with simulation volume sides of 0.8 in units where 
G = h = c = 1. d>(r) was set to 




Figure 1: A plot of the log^g of ^^e error in calculation of the elliptic scalar field equation 
for on a 75^ grid for d = 1.0. 
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(where d is a parameter) on the initial and final time surfaces. Thus, we have identical 
values of at opposite points on the initial and final surfaces. Note here that the outer 
boundary at r = r^ax = 0.8 is fairly close to r = cxd in terms of the asymptotic fields. 
This means that the action is simple to calculate at the r^ax boundary. 

We performed integrations on grids with 25^, 35^, 50^, 75^ and 100^ gridpoints. We 
also performed a few integrations on a 200^ grid, however, these calculations were very 
time consuming. 

To claim a solution, we required that the elhptical scalar field equation be zero to 
one part in one hundred thousand. We also checked that the constraint equations were 
satisfied. In general, the elhptical equation converged to 10~^ up to some value of the 
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Figure 2: A plot of the log^g of the error in calculation of the constraint equation for a' 
on a 75^ grid for d = 1.0. 



«0 



^ 40 
60 



1 1 r- 







Figure 3: A plot of the log^^Q of the error in calculation of the constraint equation for b' 
on a 75^ grid for d = 1.0. 




Figure 4: A plot of the log^^Q of the error in calculation of the constraint equation for a 
on a 75^ grid for d = 1.0. 



10 




Figure 5: The solution for field on a 75^ grid. Notice the time symmetry of the solutions 
here and in figures (^ and (|^), this indicates that our numerical methods are obeying 
the symmetry of the equations. 

parameter d. Smaller grids tended not to converge for smaller values of d than larger 
grids. Because the constraint equations are just ordinary differential equations integrated 
forwards or backwards in the r direction on the grid, they had roughly constant errors. 
In figures ([T|), (^, (j^) and (p, we plot the log^g of ^^e errors for the solutions of equations 
(p!3|), (p!4D, (|15D and (0) for a solution on a 75^ mesh with d = 1.0. Note that the error for 
(f) is everywhere less than 10~^. The errors for the constraints on a' and b' are negligible 
(parts of the plotted surface that are missing correspond to errors of to within machine 
accuracy). And the error for the constraint on d is less than 10~^'^. Significantly, the 
error for the constraint equation for d is largest where spurious oscillations occur in the 
solution (see below). 

In figures (^, (|]) and (|^ the solution for d = 1.0 is shown on a 75^ grid. Notice the 
time symmetry in (j) as well as e° and e''. This is evidence that our relaxation method 
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Figure 6: The solution for metric field e" on a 75^ grid 
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Figure 7: The solution for metric field e* on a 75^ grid 
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Figure 8: The solution for metric field e° on a 35^ grid. This is a side view. The 
amplitude of the spurious oscillations on this coarse grid decrease as the grid resolution 
increases in figures (|^) and (p!0D 



preserves the time symmetry of the equations. 

In all numerical solutions, we find spurious oscillations near, but not at the boundary, 
as presented in figures @, (|^) and ([T0|). The oscillations occur only in the a metric 



component. Notice from the figures that the amplitude of the oscillations decreases with 
finer meshes, with oscillations on the 200^ grid substantially reduced with respect to 
the coarser grids. Interestingly, the spurious oscillations do not appreciably affect the 
value of the action. This is because the oscillations are integrated over in calculating the 
action. A feature that does affect the action calculation can be seen in the depth of the 
dip near r = 0. This feature ranges in value from 0.9965 for the 35^ solution, to 0.9947 
for the 200^ solution. In figure (^TJ), we see how this influences the value of the action. 
In figure (|T2|) , we plot log^g of ^^e action as a function of log^^Q of the parameter d. The 
resulting power law with slope 2 is quite striking. In this figure, we have actually plotted 
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Figure 9: The solution for metric field e"" on a 100^ grid. 
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Figure 10: The solution for metric field e"" on a 200^ grid. 
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Figure 11: In this plot, we show the dependence of the value of the action on the 
resolution of the numerical integration. 
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Figure 12: Here, we plot log^o of the absolute value of the action as a function of log^Q 
of parameter d for gridsizes 35^, 50^ and 100^. Thus, / ~ —d?. This behavior shows the 
expected increase of persistence amplitude with binding energy of the system. 
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the absolute value of the action, which means that the action / ~ —(P. The decrease 
of the action with increasing rf is a natural result of the increased binding energy of the 
system. As we enter the nonlinear regime, we would expect that the power might change, 
and, indeed, we do begin to see a bend for larger values of d in the plot. However, the 
bend decreases with increasing grid resolution. For various grid sizes L, we find the slope 
n to be: 



(25) 



The straightening of the bend in the action is due to the better resolution of the 
dip in the metric functions on higher resolution meshes. The tendency for the slope 
to approach 2 for high resolution meshes, may indicate that the actual solution has a 
slope of 2 exactly, even well into the nonlinear regime. However, we have only shown the 
tendency, and cannot prove this assertion. 

5 Discussion and Conclusions 

In this paper, we have studied amplitudes for the spherically symmetric system of a 
scalar field coupled to gravity. We have looked at amplitudes for the field to persist in 
its current state. We have taken the field to have support primarily near the symmetry 
axis, with amplitude parametrised by the parameter d. 

What we can say from our study, is that the Euclidean action is well-behaved, negative 
and obeys obvious power law behaviour, with the action / ~ —d^. There are indications 
that the power behavior remains the same into the nonlinear regime. The decrease of 
the action with increasing d is a result of the increasing binding energy of the system as 
the scalar field gets large. This result is also familiar from the understanding of black 
hole evaporation: massive black holes radiate more than less massive black holes. The 
persistence amphtudes that we calculate show the opposite side of the coin: a massive 
scalar field is more likely to retain its integrity than a less massive scalar field. 

The solutions that we obtain are all time symmetric and largely slowly varying (for 
values of d which we have been able to probe) outside of regions near the initial and final 
times, where they vary rapidly. 

From the numerical solutions that we have obtained, it is clear that more work needs 
to be done exploring the strong-field Euclidean gravity equations. The obvious direction 
to go numerically is to develop a mesh-refinement code. A mesh-refinement code would 
increase the efficiency of the calculation by focusing on regions of the solution that are 
varying quickly, thus decreasing the computational time needed for a given resolution, 
and increasing the resolution of the solution. 

This conclusion from our work on fixed meshes is in accord with what we know of 
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critical behavior from work on the Choptuik problem. Adaptive meshes with an effective 
resolution of billions of gridpoints were necessary to see the full nonlinear critical behavior 
in the Lorentzian case. 

As well as refining the numerical approach to this problem, we plan to investigate 
more different kinds of solution. In particular, solutions which are not time symmetric, 
in order to investigate amplitudes for evaporative and condensive processes. 
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